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Abstract 

In this paper we prove the uniform boundary Harnack principle in general open sets for 
harmonic functions with respect to a large class of rotationally symmetric purely discontinuous 
Levy processes. 
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1 Introduction 

The boundary Harnack principle for classical harmonic functions is a very deep result in potential 
theory and has many important applications in probability theory and analysis. 

In the late nineties Bogdan [3] established the boundary Harnack principle for harmonic func- 
tions of rotationally symmetric a-stable processes, a € (0,2), in Lipschitz domains. This was the 
first time that the boundary Harnack principle was established for harmonic functions with respect 
to non-local operators (or, equivalently, discontinuous Markov processes). Since then the result has 
been generalized in various directions. In [T7] Song and Wu extended the boundary Harnack prin- 
ciple to harmonic functions with respect to rotationally symmetric stable processes in re-fat open 
sets, with the constant depending on the local geometry near the boundary. The definitive result 
in the case of rotationally symmetric stable processes was obtained in [4] by Bogdan, Kulczycki 
and Kwasnicki who established the boundary Harnack principle in arbitrary opens sets with the 
constant not depending on the open set itself. This type of result is known as the uniform boundary 
Harnack principle. Note that the uniform boundary Harnack principle is not true for Brownian 
motion. 

In another direction, the boundary Harnack principle has been generalized to different classes 
of discontinuous processes. In [8] the boundary Harnack principle was established for harmonic 
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functions with respect to a wide class of purely discontinuous subordinate Brownian motions in 
re-fat open sets, with an extension obtained in [ID]. In [IT] (see also, [5J[D]) the boundary Harnack 
inequality was established for harmonic functions of subordinate Brownian motions with Gaussian 
components. 

The purpose of this paper is to generalize the main results from [U [8j [ID] and prove the uniform 
boundary Harnack principle for harmonic functions with respect to a large class of rotationally 
symmetric purely discontinuous Levy processes in arbitrary open sets. The class of processes treated 
in this paper is larger than the class of processes treated in [8j [ID] . The processes considered in this 
paper need not be subordinate Brownian motions. Even when restricted to subordinate Brownian 
motions, the assumptions on the subordinate Brownian motions in this paper are slightly weaker 
than those in [51 HP]. 

To be more precise, let S = (St : t > 0) be a subordinator with Laplace exponent <\>. We assume 
that is a complete Bernstein function satisfying the following upper and lower scaling conditions 
(see [2IJ): 

(H): There exist constants #i,#2 £ (0, 1), a\-,ai > and Ro > such that 

(LSC) (j)(Xr) > ai \ Sl (t>(r), \>l,r>l/B% 

(USC) </>(\r) < a 2 X 52 (j)(r), A>l,r>l/i$. 

Note that it follows from (USC) that 4> has no drift. 

Let W = (Wt ■ t > 0) be a Brownian motion in M. d , d > 1, independent of the subordinator S. 
The subordinate Brownian motion Y = (Yj : t > 0) is defined by Yt := Ws t - The Levy measure of 
the process Y has a density given by J(x) = j{\x\) where 

/>oo 

j(r) := / (47rt)- d / 2 e- r2/(4 *V(t)^, r > (1.1) 



J o 

and fj,(t) is the Levy density of S. Note that the function r i-> j(r) is continuous and decreasing on 
(0,oo). 

We will assume that X is a purely discontinuous rotationally symmetric Levy process with Levy 
exponent v &(^). Because of rotational symmetry, the function ^ depends on |£| only, and by a slight 
abuse of notation we write ^(^) = We further assume that the Levy measure of X has a 

density Jx- Then 



e i£-(Xt-X ) 



e -**(ie]) ; for every x € M d and ^ e 



with 

*(|£|)= / (l-co S (Z-y))J x (y)dy. (1.2) 



We assume that Jx is continuous on M. \ {0} and that there is a constant 7 > 1 such that 

7 _1 j'(|y|) < Mv) < 7i(bl), for all y e R d . (1.3) 
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Clearly (|1.3p implies that 

7~ 1 <Ml£| 2 )<*(l£l)<7<Ml£| 2 ), foralUelR d . (1.4) 

For a Greenian open set D <zM d , we will use Kry to denote the Poisson kernel of X in D x D° (see 
(|4.6p below). The goal of this paper is to establish the following result: 

Theorem 1.1 Let X be a purely discontinuous rotationally symmetric Levy process with a contin- 
uous Levy density Jx satisfying (jl.3p where the complete Bernstein function 4> satisfies (H) . There 
exists a constant c = c(4>, 7, d) > such that 

(i) For every zq G R , every open set D C M. d , every r G (0, 1) and for any nonnegative functions 
u,v in M. d which are regular harmonic in D n B(zo, r) with respect to X and vanish in D c n 
B(zq, r), we have 

u(x) < u{y) 
v{x) ~ v{y) 

for allx,y G D n B(z ,r/2). 

(ii) For every zq G Mr, every Greenian open set D C M. d , every r G (0, 1), we have 

K D (x 1 ,yi)K D (x 2 ,y2) < cK D (x 1 ,y 2 )K D (x 2 ,yi) 
for all x±, X2 G D Pi B(zq,t/2) and all y±, y 2 G D° n -B(zo, ^) c - 

The proof of the above theorem uses some results developed in [lUj and several ideas from [3]. 
In the next section we recall some necessary definitions and results from [10] . In Section 3 we prove 
several results about one-dimensional symmetric Levy processes that will be needed in the proof of 
Theorem ll.il In Section 4, we present some estimates on the Poisson kernel Kd that are essential 
for the proof of Theorem II .11 The proof of Theorem II .11 is given in Section 5, where we also give an 
approximate factorization of the Poisson kernel, see Corollary 15.61 In the last section we relate the 
assumption (H) with the class OR of O-regularly varying functions and sketch the construction of 
an example of a complete Bernstein function which satisfies (H) but not the assumptions in [10j . 

At the meeting "Foundations of Stochastic Analysis" held in Banff from September 18 to 23, 
2011, M. Kwasnicki announced that, in a forthcoming joint paper with K. Bogdan and T. Kumagai, 
they have obtained a version of the boundary Harnack principle for Hunt processes in metric 
measure spaces under rather general conditions. 

In this paper we always assume d > 1. We use the following convention: The value of the 
constant C will remain the same throughout this paper, while c, ci,c 2 ,--- stand for constants 
whose values are unimportant and which may change from location to location. The dependence 
of the lower case constants on the dimension d will not be mentioned explicitly. The labeling of 
the constants ci,c 2 , ■ ■ ■ starts anew in the proof of each result. The notation f(t) X g(t), t — > 
(respectively f(t) X g(t), t — > 00) means that the quotient f(t)/g(t) stays bounded between two 
positive constants as t — > (respectively t — > 00). 
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2 Preliminaries 



Suppose that S = (St : t > 0) is a subordinator with Laplace exponent (f>, that is, S is a nonnegative 
Levy process with Sq = and 



e"^ (A) , Vt,A>0. 



E 

The function <t> can be written in the form 

-At 



0(A) = b\+ / (1 - e- M )fx(dt) (2.1) 
jo 

where 6 > and // is a measure on (0, oo) satisfying / °°(1 At)/j,(dt) < oo. b is called the drift of the 
subordinator and [i the Levy measure of the subordinator. The function is a Bernstein function, 
i.e., it is C°° and (-l) n D n <f> > for all n > 0. 

Note that, by using (|2. 1|) and the elementary inequality 1 — e~ ty < t(l — e~ y ) valid for alH > 1 
and all y > 0, we see that the Bernstein function <f> satisfies 

0(f A) < \<f>(t) for all A > 1, t > 0. (2.2) 

In this paper we will always assume that <j> is a complete Bernstein function, that is, the Levy 
measure fi of S has a completely monotone density /i(t), i.e., (—l) n D n [i > for every non-negative 
integer n. For basic results on complete Bernstein functions, we refer our readers to [16]. It follows 
from [101 Lemma 2.1] that there exists c > 1 such that 

H(t) < cn(t + 1), t> 1. (2.3) 

The next result will be used to obtain the asymptotic behavior of fx(t) near the origin. 

Proposition 2.1 ( [21|, Theorem 7]) Suppose that w is a completely monotone function given by 
w(X) = J °° e _At f(t) dt, where f is a strictly positive decreasing function. Then 

f(t) < (l-e- 1 ) -1 *- 1 ™^- 1 ), t>0. 

If furthermore, there exist 5 € (0, 1) and a, to > such that 

w(\t) < a\- 5 w{t), \>l,t>l/t , (2.4) 

then there exists c = c(f,a,to,5) > such that 

f(t) > ct^wit- 1 ), t < to. 

From now on we will always assume that the Laplace exponent <f> of S is a complete Bernstein 
function satisfying (H). 
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Theorem 2.2 For every M > 0, there exists c = c(M, eft) > 1 such that the Levy density [i of S 
satisfies 

c -1 *" 1 ^" 1 ) < fi(t) < ct~V(* _1 ) and c' 1 ^ 1 ) < fi(t, oo) < c^i -1 ), Vi < M (2.5) 
where //(i, oo) = f 00 /x(s) (is is i/te faiZ o/ i/ie Levy measure fi. 

Proof. Let w(X) := X~ 1 (j)(X) = J °° e _A */i(t, oo) (it. The upper scaling condition (USC) implies that 
w satisfies (12. 4j) with 5 = 1 — 62 and to = Rq. Hence by Proposition 12.11 there exists a constant 
ci > 1 such that 

c^t^wit' 1 ) < fi(t, 00) < cit^wO" 1 ) , t < Rl , 
which immediately implies 

qV(r ! ) < fi{t, 00) < ci(j){t~ l ) , t<i$. (2.6) 



We proceed to prove the first inequality. Since fi(t/2,co) > f t ^ 2 fi(s)ds > (t/2)/j,(t) by (|2.2p 
and for all t € (0, 

< 2i-V(*Aoo) < 2ci^ 1 </>((V2) _1 ) < 4cit _ V(* _1 )- 
Using (LSC) we get that for every A > 1 

</»(s- 1 ) = 0(A(As)- 1 )>a 1 A <5l 0((As)- 1 ), s<^. (2.7) 

Fix Ai := 2 1 / 5l ((cfa7 1 ) V I) 1 / 51 > 1. Then, by §2J3) and (HU), for s < A 1)/A X , 

/x(Ais,oo) < ci0((Ais) _1 ) < cia^ 1 \±~ Sl ^(s _1 ) < c^a^ 1 Ai _5l /i( s ) 00 ) < 2~' u ( s ' 00 ) 
by our choice of Ai. Further, 

f XlS 1 1 

(Ai — l)sfi(s) > J fi(t) dt = fi(s, 00) — /^(Ais, 00) > fi(s, 00) — 2' u ( s ' 00 ) = 2^^ S ' 00 ^ ' 



This implies that for all t < (Rq A 1)/Ai 

The case (i?§Al)/Ai < i < M is clear since the functions we consider are all positive and continuous 
on (0, 00). The proof is now complete. □ 

A consequence of (|2.5p and (USC) is that for any K > there exists c = c{K) > 1 such that 

< cn(2t), te(0,K). (2.8) 



Suppose that W = (Wt : t > 0) is a Brownian motion in R with 



E 



V£ e R a ,t> 0, 



and that W is independent of S. The process Y" = (Y t : t > 0) defined by Y t = Wg t is called 
a subordinate Brownian motion. It is a rotationally symmetric Levy process with characteristic 
exponent 3>y(£) = <K|£| 2 ), £ ^ Recall that the Levy measure of 1" has a density J(x) = j(\x\) 
with j given by (jl.lj) and that r i— > j(r) is continuous and decreasing on (0, oo). 
The following theorem establishes the asymptotic behavior of j near the origin. 



Theorem 2.3 It holds that 

j(M) ~ r nPbr^ (2-9) 



Proof. By (12T21 . 



< 0<u<t;. (2.10) 



D U 

To obtain the upper bound in ()2.9p we write 

/>r 2 />oo 

i(r)= / (4^)- d / 2 e- r2 /( 4i V(t)dt+ / (4irty d / 2 e- r2 ^ 4t ^(t)dt:= J! + J 2 . 

J0 Jr 2 

For r < 1, by using (|2.5p in the first inequality and (|2.10|) in the second, we have 

Jl<oy I" Unty^e-^h-Uit-^dtKci f (47rf)- d / 2 e^V(4t) t -i t -i^l2^ 
Jo Jo r ~ 

POD 

< c 2 r 2 <f>(r~ 2 ) / t-ip-^m dt = czr 2 <i>{r- 2 )r- d - 2 = c 3 r' d ^(r- 2 ) . 
Jo 

Next, 

/OO l-CO / J />oo \ 

t- d l 2 ^{t) dt = c 4 J I- J s-^ 2 - 1 dsj nit) dt 

fi(t) s"^ 2 " 1 (is < c 5 ^(r 2 , oo) / s^ 2 " 1 ds < cer^r" 2 ) 



= c 5 

where the last inequality follows from Theorem 12.21 The last two displays show that j(r) < 
C7r~ d (f>(r~ 2 ), for r small. To prove the converse inequality, we also use Theorem 12.21 and get that 
for r < 1, 

j(r)> l\^s)- d / 2 e- r2 l^tJL{s)ds = {±K)- d / 2 (tr 2 )- d l 2 e- r2 l^ tr2 ^{r 2 t)r 2 dt 



>c 8 r 2 - d C t- d l 2 e- l '^ l i(r 2 t)dt>^r 2 - d / r^e-V^r^-l^.^j-l) dt 



> c w r- d / t-f/s-i e -i/(*)^( r -") dt = Cll 



r "cifr 2 ) 
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where the last inequality follows because r _2 t _1 > r~ 2 and <ft is increasing. □ 
Using (|2.3p and (|2.8|) . we can easily show (see [TUl Proposition 3.5] or [131 Lemma 4.2]) that 

(1) For any M > 0, there exists c = c(M, (ft) > such that 

j(r) < cj(2r), Vr G (0, M) . (2.11) 

(2) There exists c = c((ft) > such that 

j(r) < cj(r + 1), Vr > 1. (2.12) 

3 Some results on symmetric Levy process in M. 

In this section we assume that d = 1 and denote the process X by Z. That is, (Z^F^) is a purely 
discontinuous symmetric Levy process in R such that 



e te-(z t -z ) 



e -*« r (|e|) ; for every x € R and 6> € 



We assume that (|l-4[) holds with a complete Bernstein function (ft satisfying (H), that is, r y~ l (ft{6 2 ) < 
^(1^1) < 70(0(^ 2 ) for an £ Rj but we do not assume the assumption fjl .3 j) concerning the Levy 
measure of Z. As a consequence of (H), (|1.4p and \15\ Proposition 28.1] we know that for any 
t > 0, Zt has a density pt(x, y) = pt(y — x) which is smooth. 

Let x ( K > respectively) be the Laplace exponent of the ladder height process of Z (Y, respec- 
tively). It follows from [3 Corollary 9.7] that 

It follows immediately from these two equations and (jl.4p that 7 _1 / 2 k(A) < x(A) < 7 1 / 2 k(A), i.e., 
that x is comparable to k. From (H) and \10\ Propsoition 3.7] or [12, Proposition 2.1] we conclude 
that the ladder height process of Y has no drift and is not compound Poisson, thus the ladder 
height process of Z has no drift and is not compound Poisson. Thus the process Z does not creep 
upwards. Since Z is symmetric, we know that Z also does not creep downwards. Thus if, for any 
a € R, we define 

T a = inf{t > : Z t < a}, a a = inf{t > : Z t < a}, 

then we have 

P*(r = o«) = 1, x>a. (3.2) 

Let be the process Z killed upon exiting (0, oo). Since Z has a smooth density, we can 

easily show that Z^ has a density p^>°°\t, x, y). Let G^'°°\x,y) = J °° p^'°°\t, x,y) dt be the 
Green function of Z^°°\ If we use V to denote the potential measure of the ladder height process 
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of Z, then using the symmetry of Z and pQ Theorem 20, page 176] we have that for any x € (0, oo) 
and any nonnegative function / on (0, oo) 

f OO PX 

f(y)G(°> °\x,y)dy= V(dy) V {dz)f(x + y-z). (3.3) 
o Jo Jo 

In the following, we will also use V to denote the renewal function of the ladder height process of 
Z: V(t) := V({0,t)). For any r > 0, let G (0 ' r) be the Green function of Z in (0,r). Then we have 
the following result. 

Proposition 3.1 For all r > and all x £ (0, r) 

G (0 ' r) (x, y) dy < 2V(r) (V(x) A V(r - x)) . 



Proof. Since 



f G(°< r \x,y)dy < f G^°°\x,y)dy, 
Jo Jo 



'0 Jo 

we can apply (|3.3p with / being the indicator function of (0, r) to immediately get the conclusion 
of the proposition. □ 

The following result will play an important role in this paper. 

Proposition 3.2 There exists a constant c = 0(7) > 1 such that for all r > 

r- 1 1 - < V(r) < c ' 



Proof. The proof is a simple modification of [12, Theorem 4.4]. By [101 Proposition 3.7] (or 
[121 Proposition 2.1]) we have that c^ 1 a/0(A) < k(A) < ciy^(A) for a constant ci > 1. Hence 
C 2~V<?K A ) ^ < c 2 y/4>{\), c 2 > 1, implying that 

C3 1 r—, — 7yr < CV{r) < c 3 ' 



r^/(/>(r 2 ) r-y/ </>(r 2 ) 

(where >CV(r) denotes the Laplace transform of the function V). The claim now follows by repeating 
the second part of the proof of |12[ Theorem 4.4]. □ 



4 Poisson Kernel Estimates 

Recall that Y is a subordinate Brownian motion in M. d with Levy exponent <^(|£ | 2 ), X is a purely 
discontinuous rotationally symmetric Levy process in M. d with Levy exponent = ^(1^1) and 

Levy density Jx, i.e., 



E 3 



-tm\) 



for every x € 



l d and 
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and ^(|^|) = f™ d (l — cos(£-y))Jx(y)dy. Recall that we assume that (|1.3p holds. As a consequence of 
(H), (jl.4p and [151 Proposition 28.1] we know that for any t > 0, X t has a density pt(x, y) = pt(y—x) 
which is smooth. 

The infinitesimal generator L of X is given by 

L/(z) = / {f(x + y)- f{x) - y ■ V/(x)l {|l/ |< 1} ) Jx(y)dy (4.1) 

for / € Cf(R d ). Moreover, for every / € C£(R d ), /(X t ) - /(X ) - f* Lf(X a ) ds is a P x -martingale 
for every x € 

First we record several inequalities that will be needed in the remainder of the paper. 
Lemma 4.1 There exists a constant c = c((f>) > such that 

<P(r- 2 ) 1/2 dr < cA-VCA 2 ) 172 , VA > 1/R , (4.2) 



o 



f-A- 1 rR 

A 2 / rct){r- 2 )dr+ r"V(r~ 2 )dr < c<^(A 2 ), VA > l/i? , (4.3) 



JA 

and 

r(j)(r- 2 ) 1/2 dr+ r~ x <\>{r- 2 fl 2 dr < cc/>(A 2 ) 1/2 , VA > l/i?o- (4.4) 

/o A- 1 



v2 



Proof. Assume A > 1/Rq. By (USC), 0(r~ 2 ) < c ir - 2S2 \- 2 ^<f>(\ 2 ) for r < A -1 . On the other 
hand, by (LSC), 0(r~ 2 ) < c 2 r- 2Sl \- 2S ^{\ 2 ) for A" 1 <r< R . Thus 

" 1 0(r- 2 ) 1 / 2 d r < c ;/ 2 ^(A 2 ) 1 /2 A -«5 2 1 r -^ dr < C3 A-i0 (A 2 ) i/2^_ ) 

7o 1-02 

/•A" 1 pRo 
A 2 / rcf)(r~ 2 )dr + I r~ l (f)(r~ 2 )dr 
Jo Jx- 1 

< c 4 0(A 2 ) ^A 2 - 252 £ 1 r 1 " 25 ^ + r" 1 " 2 * 1 ^ < c 5 0(A 2 ) f ^ 1 - , + 1 



2(1 -<5 2 ) 2*i 



and 



/•A- 1 /-iio 
A 2 / r(\>(r- 2 f t2 dr + / r -1 0(r~ 2 ) 1/2 dr 

io JA- 



-i 



< cg^A 2 ) 1 / 2 L 2 - 52 jf * 1 r 1 "^ + A" 51 1^ r'^dr] < c^X 2 ) 1 ' 2 (-^- + ±) 



□ 
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Lemma 4.2 There exists a constant c = c(</>, 7) > such that for every f G C%{R d ) withO< f < 1, 
L/ r (x) < c0(r- 2 ) [2 + isup^|(5 2 /9y^y fc )/(y)| ] + b , for every x e R d , r < R 



3,k 



where f r (y) := f(y/r) and b := 2 f\ z \ >Ro Jx{z)dz < 00. 

Proof. Let L x = su Py \(d 2 /d yj dy k )f(y)\. Then \f(z + y) - f(z) - y ■ Vf(z)\ < ^L^ 2 . For 
r G (0,-Ro]) let f r (y) = f(y/r). Then the following estimate is valid: 

L\ \y\ 2 

1/rO +y) - fr(z) - y ■ V/ r (z)l { |^|< r} | < — 2"l{|y|<r} + 2 ■ l{|y|> r } • 

Now, by using (H), ([2JJ]) and we get 

\Lfr(z)\ < / \fr(z + y)~ fr(z)-y-Vf r (z)ln y \< r y\Jx{y)dy 

Jm d 

- ~T I 1 {\y\<rY^T J x{y)dy + 2 / l {r <\ y \< Ro} J x (y)dy + 2 [ ln y \> Ro xJ x {y)dy 

1 JR d r JR d JR d 



< -> L 



-^T I 1 {\y\<rY^\rj{\y\) d V + 2 1 I 1 {r<\y\<R }j(\y\)dy + 2 f l { \ y \> Ro} Jx{y)dy 

z JR d r JR d JR d 



where the constant c is independent of r £ (0,i?o]- n 



Jx(y)dy, 
'{\y\>Ro} 



For any open set D, we use T£> to denote the first exit time of D, i.e., td = inf{i > : Xt ^ I?}. 
Using Lemma [4.2l the proof of the next result is the same as those of |10} Lemmas 4.1 and 4.2]. 
Thus we skip the proof. 

Lemma 4.3 There exists a constant c = 0(^,7) > such that for every r G (0,1], and every 
x G M d , 



The idea of the proof of the following proposition comes from |19j . 

Lemma 4.4 There exists c = 0(7) > such that for any r G (0, 00) and xq G M. d , 

Ex[TB(x ,r)] < c (0(r- 2 )^((r - fx - xqI)- 2 ))- 1 ^ s g r ). 

Proof. Without loss of generality, we may assume that xq = 0. We fix x ^ and put Z t — 
Then, using the fact that ^ is a radial function, Zt is a Levy process on R with 

E[e^} = E(e t6 W Xt ) = e^wO = e -**W Q G R. 
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Clearly, 7~ 1 (/>(0 2 ) < *&(0) < 7</>(# 2 ). Thus Zt is of the type of one-dimensional symmetric Levy 
processes studied in Section [3J 

It is easy to see that, if X t G .6(0, r), then \Z t \ < r, hence E a ,[r B ( 0r )] < Euiff], where f = 
inf{t > : \Zt\ > r}. Thus, applying Proposition 13-H we obtain E x [T£(o,r)] < 2V(2r)V(r — \x\). 
Now, by Proposition 13.21 and (H), we have proved the lemma. □ 

We now recall the definition of harmonic functions with respect to X. 

Definition 4.5 Let D be an open subset of~R d . A function u defined on M. d is said to be 

(1) harmonic in D with respect to X if 

E x [\u(X TB )\] < oo and u(x) = E x [u{X TB )\ , x £ B, 
for every open set B whose closure is a compact subset of D; 

(2) regular harmonic in D with respect to X if it is harmonic in D with respect to X and for each 

x € D, 

u(x)=E x [u(X TD )]. 

Since our X satisfies [5, (1.6), (UJS)], by [5, Theorem 1.4] and using the standard chain 
argument one have the following form of Harnack inequality. 

Theorem 4.6 For every a S (0,1), there exists c = c(o, 0, 7) > such that for every r € (0,1), 
xq G IR d , and any function u which is nonnegative on M. d and harmonic with respect to X in 
B(xq, r), we have 

u[x) < cu(y), for all x,y E B(xo,ar) . 

Given an open set D C M. d , we define XP(lo) = X t (oj) if t < td(lu) and XP(lo) = d if t > td(w), 
where d is a cemetery state. A subset D of M is said to be Greenian (for X) if X D is transient. 
When d > 3, any non-empty open set D C M. d is Greenian. An open set D C M. d is Greenian if 
and only if D c is non-polar for X (or equivalently, has positive capacity with respect to X). In 
particular, every bounded open set is Greenian. 

Since X has a smooth density, using the strong Markov property, it is standard to show that 
for every Greeninan open set D, X^ has a density pr>(t,x,y). For any Greeninan open set D in 
M. d let Go(x,y) = J* °° pn{t,x,y) be the Green function of X D . Using the Levy system for X, we 
know that for every Greeninan open subset D and every / > and x € D, 

E x [f(X TD ); X TD _ ± X TD ] = [_ [ G D (x,z)J x (z - y)dzf{y)dy. (4.5) 

Jd c Jd 

We define the Poisson kernel 

K D (x,y) := G D (x,z)J x (z -y)dz, (x,y)eDxD c . (4.6) 
Jd 
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Thus (|4.5p can be simply written as 

E x [f(X TD ); X TD _ + X TD ] = [_ K D (x,y)f(y)dy. 

Jd c 

Using continuity of Jx, one can easily check that Kd(%, •) is continuous on D c for every x G D. 

Proposition 4.7 There exists c\ = ci(4>, 7) > and C2 = c%((f), 7) > such that for every r G (0, 1] 
and xo £ 

K B (xo,r){v,y) < cij(\y-x \-r)(<j>(r- 2 )<f>((r-\x-xo\)- 2 )y 1/2 (4.7) 
< cudy-xol-r)^- 2 )- 1 (4.8) 



for all (x,y) G B{xQ,r) x B(xo,r) and 

K B (x ,r)(.x ,y) > c 2 j(\y - xoD'K'" -2 ) -1 , for all y G B(x ,rf. (4.9) 

Proof. Using (TOj) and (pmJ)Hf2TT2j) . the proof of P~Tj) and (jM]) is exactly the same as that of 
[10\ Proposition 4.10] (using (H)), while (|4.8p follows from (|4.7p and the fact that is increasing. 
□ 



Using Theorem 14.61 and the continuity of K B / XQir \(x, •) on B(xo,r) for every x £ D, the proof 
of the next result is the same as that of [10\ Proposition 1.4.11]. So we omit the proof. 

Proposition 4.8 For every a G (0, 1) ; there exists c = c(4>, 7, a) > such that for every r G (0, 1], 
xq G M. d and xi,x 2 G B(xq, ar), 



K B(x ,r){xi, y) < cK B{x0!r) (x 2 ,y), y£B(x ,r) . 

Proposition 4.9 For every a G (0, 1), there exists c = c(</>,7, a) > such that for every r G (0, 1] 
and xq £~R d , 

K B(x(hr) {x,y) < c r ^ ^— ^ j 

for all x G B(xo, ar) and all y such that r < \xq — y\ < 2r. 
Proof. By Proposition 14.81 

K B{x ,r){x,y) < % I K Bix0tr) (w,y)dw 
for some constant c\ = ci((p,j,a) > 0. Thus from Lemma [4.41 (|4.6p and Theorem 12.31 we have that 
K B (x (U r)(x,y) < % / / G B(x0ir) {w, z)J x (z - y)dzdw 

r JB(x ,r) JB(x ,r) 



rd JB(x ,r) 



^z[T B ( X(ur) ]Jx(z - y)dz 



< 02 f <K\*-v\~ 2 ) \ z - v \-d dz 

~ r d {4>(r- 2 )Yl 2 J B[xo , r) (0((r - \z - xqI)" 2 )) 1 ^ 1 y] 
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for some constant c 2 = 02(^,7, a) > 0. Since r — \z — xq\ < \y — z\, we have 



K B (x , r )(x,y) < ,/! 2m/2 / — — dz 

r d {(f>{r 2 )y/ 2 J B (x ,r) \z-y\ d 



{<j>{\z-y\- 2 ))V\ 

B(xo,r) 



< /" OKI* -^l" 2 )) 1/2 d - 

r d (0(r- 2 ))V2 7B(v,3r)\B(y,| W -*o|-r) I 2 _ J/| d 



3r <Xs- 2 ) 1/2 



Now, using (|4.4p in the last integral (considering the cases r < i?o/3 and 1 > r > Rq/3 separately), 
we arrive at the conclusion of the proposition. □ 



Lemma 4.10 For every a G (0, 1), there exists a positive constant c = c(</>, 7, a) > such that for 
any r G (0, 1) and any open set D with D C -6(0, r) we have 

F X (X TD G#(0,r) c ) < C 0(r" 2 ) / G D (x,y)dy, x G D n B(0,ar) . 

Jd 

Proof. The proof of the lemma is similar to that of \10\ Lemma 4.15]. Transience was used in 
the proof |1U[ Lemma 4.15] in order to derive equation \10\ (4.29)]. By noting that \10\ (4.29)] in 
the proof of [10, Lemma 4.15] follows immediately from Dynkin's formula, using our Lemma 14.21 
we can follow the rest of the proof \10\ Lemma 4.15] (which does not use transience) to get the 
conclusion of the lemma here. We omit the details. □ 



5 Uniform Boundary Harnack Principle 

In this section, we give the proof of the main result of this paper. Let A(x, a, b) := {y G M. d : a < 
\y — x\ < b}. 

Lemma 5.1 For every p G (0, 1), there exists c = c((p,^,p) > such that for every r G (0, 1), 

■\y\ r 

K B{0tS) (x,y)ds < c-—^j(\y\) Vx G B(0,pr), y G A(0,r(l +p)/2,r). 

r(l+p)/2 n r ) 

Proof. Let < p < 1 and q = (1 + p)/2. Note that the functions r 1-4 r~ d+1 and r 1— >■ 
r -1 ^^ -2 )) -1 ' 2 are decreasing, see ()2.10p . Using Proposition 14.91 we get 
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for some constants c±(p, (j>) > and ozip, <j>) > 0. Note that by (|4.2p (considering the cases \y\ — qr< 
Rq and 1 > \y\ — qr > Rq separately) and the fact that r i— > r(0(r -2 )) 1 / 2 is increasing 

r\y\ r\v\~i r 

/ m\y\-s)- 2 )) 1/2 ds = / {^(s-^ds < cM-qrmdyl-qr)- 2 )) 1 / 2 < c 3 r(0(r- 2 )) 1 /2 

Jgr JO 

for some constant C3 > 0. Thus, by (H), Theorem 12.31 and the fact that r — > j(r) is decreasing, we 
have 

y\ r r 



' K B{0)S) (x,y)ds < < c 5 ^Z2ji(r) < c 6^2)j(M)- 



□ 



From the strong Markov property, it is well known and easy to see that for every Greenian open 
sets U and D with U C D, G D (x,y) = Gu(x,y)+E x [G D (X Tu ,y)] for every (x,y) € R d xR d . Thus, 
for every Greenian open sets U and D with U C D, 

K D (x,z) = K u (x,z)+E x [K D (X TU ,z)], (x,z)eUxD c (5.1) 

and 

[ TD ] = E x [ m ] + M x [e Xtu [r D ]] x £ V. (5.2) 

Lemma 5.2 For every p £ (0, 1), i/iere exists c = c(0, 7,p) > suc/i i/iai /or every r € (0, 1), /or 
every zo € R d , [7 C B(zo,r) and for any (x,y) € (U C\ B(zo,pr)) x S(zo,r) c , 

Ku(x,y) < c ( / j(l z - zo\)Ku(z,y)dz +j(\y - z \) ) . 

<P{r ) \Ju\B(z ,(l+p)r/2) J 

Proof. Without loss of generality, we assume zq = 0. Let < p < 1, qi := (1 +p)/2 and 
<Z2 := (3 + 2p)/5. For every s G [ojr, q^r] and x G C7 n B(0,pr), by (|5.ip we have 

Ku(x,y) = E x [K u (X TunB(0s) ,y)}+ K UnB(0tS) (x,y) 

Ku(z,y)K UnB{ Q >s )(x,z)dz + ^f/oBfo^O^S/) 

(7\S(0,s) 

< / K u (z,y)K B ( 0>s - ) (x,z)dz + K B ( 0>s - ) (x,y). 

JU\B(0,s) 



Thus 



Ku(x,y)<— -/ / K B ( 0s) (x,z)K u (z,y)dzds + — -/ K B(0s) (x,y)ds 

r{q2 - qi) J qir Ju\B(a, s ) r m - qi) J qir 
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By Tonelli's theorem, we have 



10 f q2T f 

r(l-p) J qir J{ z eU;\z\> qi r} 



10 

< 



r(l~p) J(U\B(0, qi r)) 



K B (o,s){x,z)ds Kjj(z,y)dz. 



qir 



Applying Lemma 15. II to the inner integral above, we get that 

ci 1 

(1 -p) <j){r- 2 ) J(v\B{0 m r)) 

One the other hand, for any s E [qir, q 2 r], by Proposition 14.71 



I< n ± j j(\z\)K u (z,y)dz. (5.3) 



K B{0 , s) (x,y) < c 23 (\y\ - «) (0(a _ 2))1/2 _ ^ y2))1/2 ■ 

When y E A(0,r,4) we have (1 — q 2 )\y\ < \y\ — s, while when \y\ > 4 we have \y\ — s > \y\ — 1. 
Since s — \x\ < s < q 2 r, we have by the monotonicity of j, 

and 

m - s) w^W 2 m - l\)- 2 )) i/2 - C3J(M ~ l) w^r !y| - 4 

for some constant C3 > 0. Thus by applying (|2.1ip and ()2.12p . we get 

" ~ C4(1 - p) L m ~ s \o(s-^ms-w)-^ ds £ c " m) W^y (5 ' 4) 



Combining (|5.3p - ()5.4p . we conclude that 



If 1 

Ku(x,y) < cq— — 2T / j(\z\)K u (z,y)dz + c G j(\y\)— — 

4>{r 2 ) J(u\B(o, qi r)) <P{r V 



□ 



Note that, since X satisfies the hypothesis H in [18], by [181 Theorem 1], if V is a Lipschitz 
open set and U C V 

P X (X TU E dV) = and F x (X Tu E dz) = Ku(x,z)dz on V c . (5.5) 

Lemma 5.3 For every p £ (0, 1), i/iere exists c = c(<p, J,p) > such that for every r E (0, 1), for 
every zq E M. d , U C B(zo,r) and any nonnegative function u in R rf which is regular harmonic in U 
with respect to X and vanishes in U c n B(zo,r) we have 

u(x) < c [ j(\y - z \)u(y)dy, x E U n B(z ,pr). 

<P\ r ) J (U\B(z ,(l+p)r /2))US(z ,r-) c 
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Proof. Without loss of generality, we assume zq = 0. Let < p < 1 and set q = (1 +p)/2. 
Note that the part of boundary of U belonging to U c n B(zo,r) needs not be Lipschitz, but here 
u vanishes. The other part of boundary of U is a part of the boundary of the ball B(zq,v). Thus, 
since u is regular harmonic in U with respect to X and vanishes in U c H B(zQ,r), by Lemma [ 
and (15. 5ft we have 

u{x) = E x [u(X TU )] = / K u (x,y)u(y)dy 

Ju c 

< c — 1 ( / / j(\z\)Ku(z,y)dzu(y)dy + I j(\v\)u(y)dy 

9V ) \JUc JU\B(0,qr) JB{0,rY , 

Since J Uc K\j{z, y)u{y)dy = u(z) on U \ B(0, qr), by Tonelli's theorem, we have 

1 



uix) < c- 



j(\z\)u(z)dz + j(\y\)u(y)dy . 

•) JB(0.r) c / 



r ) \Ju\B(0,qr) JB{0,r) 



□ 



Lemma 5.4 There exists C = C(</>, 7) > 1 such that for every r G (0,1), /or every zq G 
f7 C B(z ,r) and for any (x, y) £ U D B(z ,r/2) x (-B(z , r) c n TJ C ), 

j(\z - z \)Ku(z,y)dz + j(\y - z \) 

U\B(z ,r/2) I 



\Jv 

< Ku(x,y) <CE x [ru) I [ 



j(\z - z \)Ku(z,y)dz + j(\y - z Q \) . 

U\B(z 0} r/2) J 



Proof. Without loss of generality, we assume zq = 0. Fix r G (0, 1) and let B := B(0,r), U\ : 
±r), U 2 :=Un B(0, |r) and U 3 := U D B(0, § 



Ur\B(0±r), U 2 :=UnB(0,h) and U 3 := UnB{0, §r). Let x G Ur\B(0,r/2), y G B(0,r) c n[/ c . 



By HSU), 

= / Ku(z,y)P x (X T £dz)+ K u (z,y)K U2 (x,z)dz + K U2 (x,y) 

Ju 3 \u 2 Juxuz 



Ku{z,y)f x (X TU edz) + / Ku(z,y) / G{7 2 (x,u;)j(|2; - w\)dwdz 
u 3 \u 2 Ju\u 3 Ju 2 



+ / G U2 (x,w)j(\y - w\)dw =: I + 11 + III. 

JUo 



HJ 2 

From Lemmas 14.101 and 15.21 we see that there exist Cj = q(^>, 7) > 0, i = 1, 2, such that I is less 
than or equal to 



ci ( sup Ku(z,y) ) (j){r 2 )E x [t U2 ] < c 2 E x [tu 2 ] \ 
\z£U 3 J \Jv 



j(\z\)Ku(z,y)dz+j(\y\) . (5.6) 
u\u 3 J 
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Note that if 4 > \z\ > |r and \w\ < |r, then 

!, , 8. . , , 2 . , , . . 2 

— \z\ = \z\ \y\ < \z\ r < \z\ — \w\ < \z — w < m + u> < \z\ H — r < 2\z\. 

Thus, since j is monotonely decreasing, by (|2.1ip 

1 3 2 

C3 ^(kl) < i( 2 M) < id* - w l) < i(gkl) < c 3 j(\z\), if 4 > |z| > -r, \w\ < -r 

for some constant C3 > 0. If 4 < |z| and < |r, then |z| — |r < |z — io| < \z\ + |r and by (|2.12|) 

c^jd^l) < C4 1 j(|z|+-r-l) < j(\z\+-r) < j(\z-w\) < j(\z\--r) < c 4 j(|z|--r+l) < c 4 j(|z|), 
for some constant c 4 > 0. Thus there exists C5 = 05(^,7) > 1 such that 

c^E x [r U2 ] [ j(\z\)Ku(z,y)dz <II <c 5 E x [t U2 ] [ j(\z\)K v (z,y)dz (5.7) 
Ju\u 3 Ju\u 3 

and 

c^E x [T U2 }j(\y\) < III < c 5 M x [r U2 ]j(\y\) . (5.8) 

Now the upper bound follows from (|5.6[) — f|5.8j) . To prove the lower bound we can neglect /. Further, 
by using Lemmas 14.41 and 14.101 in the third line, from (15. 2p we get 

E x [ru] = E x [t U2 ]+E x \e Xtv [ru] 



< E x [t U2 ] + sup E* fa] f x {X G B(0, 2r/3) c 

Kzeu 



< E x [t U2 } + cs^r- 2 )- 1 C7<p((2r/3)- 2 )E x [T U2 } < c 8 E x [t U2 
for some constants cs > 0. In the last inequality above we have used (H). Since 



j(\z\)Ku(z,y)dz = / j(\z\)Ku(z,y)dz + / j(\z\)K u {z,y)dz 
U\Ui Ju\u 3 Ju 3 \U! 



< / j(\z\)Ku(z,y)dz +[ sup Ku(z,y) ) j(\y\)dy, 

JU\U 3 \z£U 3 J JA(0,r/2,3r/4) 

by Theorem 12.31 and Lemma 15.21 

j(\z\)Ku(z,y)dz < (1+-^—/ S -^(s- 2 )ds ) { I j(\z\)K u (z,y)dz+j(\y\) 



U\U! \ n r 2 ) Jr/2 I \JU\U. 



r" s - i <t>{s- 2 )ds \ ( f 

Jr/2 J \JU' 



Applying (|4.3p (considering the cases r < 4i? /3 and 1 > r > 4i? /3 separately) and (H), we 
obtain 



j(\z\)K u (z,y)dz<c 10 [ j(\z\)Ku(z,y)dz + j(\y\) . (5.9) 

U\Ui \Ju\u 3 J 

Combining (15,7p - (15,9p . we have proved the lower bound. □ 
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Lemma 5.5 For every zq G M. d , every open set U C B(zq, r) and for any nonnegative function u 
in M. d which is regular harmonic in U with respect to X and vanishes a.e. in U c D B(zo,r) it holds 
that 



C- 1 K x [tjj] / j(\y - z \)u(y)dy < u(x) < CE x [t v ] / j(\y - z \)u(y)dy 

JB(zo,r/2) c JB{z ,t/2Y 

for every x G U (~l B(zo,r/2) (where C is the constant from Lemma \5.1$ . 

Proof. Without loss of generality we may take zq = 0. By the argument in the proof of Lemma 
and by the assumption that u vanishes a.e. on U c n B(0, r) we have that 



u(x) = / Kjj(x,y)u(y)dy = I Ku(x,y)u(y) dy . 

JU c JB(z ,r) c 

Now the claim follows from Lemma 15.41 Indeed, by Tonelli's theorem we get 

u(x) < CE x [tu]I f j{\z\)l f K u (z,y)u(y)dy)dz+ [ j(\y\)u(y)dy 

\JU\B(0,r/2) \JB(0,r) c J J B(0,r) c y 



CE x [ Tu ]\ I j(\u\)u(z)dz+ j(\z\)u(z)dz 

IU\B(0,r/2) J_B(0,r) c 



= CE x [ru] / j{\u\)u{z)dz, 

JB(0,r/2) 

where for the last line we used that u vanishes a.e. on U c n B(0, r). The lower bound follows in the 
same way. □ 

We remark that in the statements of Lemmas I5.1H5.51 by using the assumption (jl.3p , we could 
have replaced the density j with the density Jx of the process X (with a different constant). We 
will do this in the next corollary which gives an approximate factorization of the Poisson kernel. It 
is an immediate consequence of the last two lemmas. 

Corollary 5.6 Let zq £ M. d , D C M. d be Greenian open set and denote U := D n B(zQ,r). Then 
for every r G (0, 1) and all (x,y) G (D l~l B(z ,r/2)) x (D c D B(z ,r) c ) it holds that 

C- l E x [Tu)A(y) < K D (x,y) < CE x [Tu]A(y) , (5.10) 

where 

A(y):= (j x (z-z )K u (z,y)dz + Jx(y-z ))+ J x (z- z )E z [K D (X TU ,y)\ dz . 

JU\B(z ,r/2) JB(zo,r/2) c 

Proof. Without loss of generality, we assume zq = and D n B(0, r/2) ^ 0. We first note that by 
(I5TTD and ([53]), for every (x, y) G {D n B{0, r)) x {D c n B(0, r) c ), 

K D (x, y) = Ku(x, y) + E x [K D (X TU ,y)\ . 
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The function x i— >■ E x [ifo (X Tc/ , y)] is regular harmonic in U with respect to X and vanishes a.e. in 
U c n -B(0, r). By using Lemma 15.51 for this function, and Lemma 15.41 for Kjj(x,y) we immediately 
obtain required inequalities. □ 



We can now easily prove Theorem 11.11 
Proof of Theorem 11.11 (i) This follows immediately from Lemma 15.51 with c := C 4 . 
(ii) Let xi,x 2 € D n B(z , r/2), yi,y 2 e» c n%r) c and let [/:=DnB(z ,r). Then by (foTTUD 

K D (x 1 ,y 1 )K D (x 2 ,y 2 ) < (CE Xl [T U ]A(yi)) (CE X2 [T U ]A(y 2 )) 

= (CEMAiy^iCEMAiyx)) 
< C 4 K D (xi,y 2 )K D (x 1 ,y 2 ) . 

The lower bound is proved in the same way. □ 



6 Remarks on (H) 

In this section we point out the relationship between the assumption (H) and the class OR of O- 
regularly varying functions, and sketch the construction of a complete Bernstein function <p which 
satisfies (H) but not the assumption in [10J that <j> is comparable to a regularly varying function. 
Using the idea in the construction below, one can come up with a complete Bernstein function that 
is bounded between any two regularly varying complete Bernstein functions. 

It follows from the definitions on [21 page 65 and page 68] and [21 Proposition 2.2.1] that the 
assumption (H) is equivalent to that <f> is in OR with its Matuszewska indices contained in (0, 1). 

Here is a sketch of the construction. For x £ (0, 2], define 

f(x)=x 1 / 2 . 

Then we define 

f{x) = x 1 / 3 + /(2) - 2 1 / 3 , x€(2,ai] 

for some large constant a\ > 2. The constant a\ is chosen so that for large values of x in (2, ai], 
the function / behaves like x 1 / 3 , that is f{Xx)/f{x) is close to A 1 / 3 uniformly for A € [1,2]. Then 
we define 

f{x) = x l/2 + f(a 1 ) - a 1 / 2 , ace(ai,a 2 ] 

for some large constant a 2 > a%. The constant a 2 is chosen so that for large values of x in (ai, a 2 ], 
the function / behaves like x 1 / 2 , that is f(Xx)/f(x) is close to A 1 / 2 uniformly for A € [1,3]. Then 
we define 

f(x) = x 1/3 + f(a 2 ) - 4 /3 > xe(a 2 ,as] 
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for some large constant 03 > a 2 ■ The constant 03 is chosen so that for large values of x in (a 2 , 03] , 
the function / behaves like x 1 / 3 , that is f(Xx)/f(x) is close to A 1 / 3 uniformly for A E [1,4]. We 
repeat this procedure to define this function inductively. 

The function / is an increasing function in OR with upper Matuszewska index 1/2 and lower 
Matuszewska index 1/3. 

Let a be the measure with distribution function /. Since /( Ooo )(l + t)~ 1 a(dt) < 00, a is a 
Stieltjes measure. Let 

J(Q,oo) A + 1 

be the corresponding Stieltjes function. It follows from integration by parts that 
Using our construction of / we know that 

tl/3 roo ff£\ ,-oc el/2 



(A+e) 2 s -y 2 (A+0 2 s -y 2 (A+0 

Thus it follows from [20, Lemma 6.3] that 

ciA^ 2/3 < 5 (A) < c 2 A" 1/2 , A > 2 

for some positive constants c\ < c 2 . 

Modifying the argument of the proof of the de Haan-Stadtmiiller theorem ([21 Theorem 2.10.2]) 
one can show that g is in OR with upper Matuszewska index —2/3 and lower Matuszewska index 
— 1/2. It follows from [16^ Theorem 7.3] that 4>(x) := l/f(x) is a complete Bernstein function. 
Thus (f) is a complete Bernstein function in OR with upper Matuszewska index 1/2 and lower 
Matuszewska index 1/3. 

It follows from [20\ Lemma 6.3] that g cannot be comparable with any regularly varying function 
at infinity, and therefore (ft cannot be comparable with any regularly varying function at infinity. 

Acknowledgements. We thank Nick Bingham and Charles Goldie for their help related to 
the construction of the function in the last section. 
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